Mon. Not. R. Astron. Soc. 000. ITHTOl (2009) Printed 22 July 2009 (MN WI$i style file v2.2) 



Follow-up observations at 16 and 33 GHz of extragalactic 
sources from WMAP 3-year data: I — Spectral properties" 

^ ! AMI Consortium: Matthew L. Davies, 1 ! Thomas M. O. Franzen, 1 ^ 
^ '. Rod D. Davies, 2 Richard J. Davis, 2 Farhan Feroz, 1 

Ricardo Genova-Santos, 1 ' 3 Keith J. B. Grainge, 1 ' 4 David A. Green, 1 
Michael P. Hobson, 1 Natasha Hurley- Walker, 1 Anthony N. Lasenby, 1 ' 4 
Cn| ■ Marcos Lopez-Caniego, 1 Malak Olamaie, 1 Carmen P. Padilla- Torres, 3 
Guy G. Pooley, 1 Rafael Rebolo, 3 Carmen Rodrfguez-Gonzalvez, 1 
Richard D. E. Saunders, 1 ' 4 Anna M. M. Scaife, 1 Paul F. Scott, 1 
£j ! Timothy W. Shimwell, 1 David J. Titterington, 1 Elizabeth M. Waldram, 1 
^ ! Robert A. Watson 2 and Jonathan T. L. Zwart 1 

O ii 1 Astrophysics Group, Cavendish Laboratory, 19 J. J. Thomson Avenue, Cambridge CBS 0HE 

I 2 Jodrell Bank Centre for Astrophysics, University of Manchester, Manchester M13 9PL 

| ^Institute de Astrofisica de Canarias, 38200 La Laguna, Tenerife, Canary Islands, Spain 

-4_> , 4 Kavli Institute for Cosmology Cambridge, Madingley Road, Cambridge, CB3 OH A 

a 



o 



(N 
> 

o 



Accepted ????. Received ???? 



ABSTRACT 

We present follow-up observations of 97 point sources from the Wilkinson Microwave 
Anisotropy Probe ( WMAP) 3-year data, contained within the New Extragalactic 
WMAP Point Source (NEWPS) catalogue between -4° < 5 < 60°; the sources form 
t a flux-density-limited sample complete to 1.1 Jy (~ 5cr) at 33 GHz. Our observations 

£^ , were made at 16 GHz using the Arcminute Microkelvin Imager (AMI) and at 33 GHz 

with the Very Small Array (VSA). 

94 of the sources have reliable, simultaneous - typically a few minutes apart - 
observations with both telescopes. The spectra between 13.9 and 33.75 GHz are very 
different from those of bright sources at low frequency: 44 per cent have rising spectra 
( a ii'9 5 < 0.0), where S oc v~ a , and 93 per cent have spectra with etfjj'g 5 < 0.5; the 
52 ' median spectral index is 0.04. 

For the brighter sources, the agreement between VSA and WMAP 33-GHz flux 
densities averaged over sources is very good. However, for the fainter sources, the VSA 
tends to measure lower values for the flux densities than WMAP. We suggest that the 
main cause of this effect is Eddington bias arising from variability. 



o 



1 INTRODUCTION 

Extragalactic point sources contaminate maps of the cosmic 
microwave background (CMB), such as those produced us- 
ing data fr om the Wilkinson M icrowave Anisotropy Probe 
(WMAP) (Be nnett et al.ll2003 ). particularly at frequencies 
< 60 GHz. Catalogues of contaminating point sources are, 
therefore, a natural by-product of CMB surveys and are an 

* We request that any reference to this paper cites 'AMI Consor- 
tium: Davies et al. 2009' 
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invaluable resource in the study of point sources at high ra- 
dio frequencies, the statistical properties of which remain 
relatively unknown. 

The WMAP mission has produced all-sky maps of the 
CMB in five frequency bands between 23 and 94 GHz. The 
point-source catalogue resulting from the original analysis o f 
the 3-year data contains 323 entries l|Hinshaw et alj [20071 ): 
spectral data from WMAP suggest that a large propor- 
tion of the detected sources have flat spectra, implying that 
they belong to a po pulation of compact sources (see e.g. 
iToffolatti etallll998l h 

lLopez-Caniego et alj (|2007l ) have re-analysed the 3-year 
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WMAP data to produce the New Extragalactic WMAP 
Point Source (NEWPS) catalogue of extragalactic point 
sources, which contains 369 sources detected at 5a in at 
least one of the frequency bands. The number of sources 
in the NEWPS catalogue has been maximised by using 
prior knowledge of so urce positions from low-frequency data. 
iHinshaw et al.l {20071 ) neglect the deviations, increasing with 
frequency, of the WMAP point spread function from a Gaus- 
sian; the NEWPS estima tes use the real beam s hape at ev- 
ery frequency. Moreover, lLopez-Caniego et alj {2007T ) have 
applied a method to correct flux density estimates for the 
Eddington bias, which leads to an overestim ate of the flux 
densities of faint sources (see IWand |2004 and references 
ther ein) . 

iDe Zotti et all (|2005h made some predictions for high- 
frequency radio surveys of extragalactic sources. In par- 
ticular, they claim that the dominant population at 
S30GHZ > 1 Jy are blazars in which the radio emission is due 
to a relativistically beamed jet observed end-on. However, 
spectral studies of complete samples selected at cm wave- 
leng ths are few. 

iBolton et alj {2004) made simultaneous VLA observa- 
tions at 1.4, 4.8, 22 and 43 GHz o f sources selected from the 
9C survey {Waldram et al.ll2003l) at 15 GHz which has a 
completeness limit of ta 25 mjy. IClearv et all (|2005l . 120081 ) 
studied a different sample from the 9C survey, this time 
with a completeness limit of « 20 mjy at 15 GHz, using 
the Very Small Array (VSA) source subtractor at 33 GHz. 
iMassardi et alj (120081 ) have made simultaneous observations 
with the Australia Telescope Compact Array (ATCA) at 4.8 
and 8.6 GHz of a sample of sources, complete to « 0.5 Jy at 
20 GHz, with the ATCA. As a next step, we here address 
the issue of source spectra at cm wavelengths of higher flux 
density sources. 

We have carried out observations of sources contained 
in the NEWPS catalogue using the Arcminute Microkelvin 
Imager (AMI) at 16 GHz and the Very Small Array (VSA) 
at 33 GHz. For the purpose of obtaining spectra, sources 
were observed simultaneously (typically within a few min- 
utes) with the AMI and VSA. Finally, we compared WMAP 
flux densities with AMI/ VSA flux densities, investigating 
the effects of source variability. 

The layout of this paper is as follows. The remain- 
der of this section describes the NEWPS Catalogue and 
our source sample. In Section [5] we describe the AMI and 
VSA telescopes. The observations and data reduction are 
explained in Section|3] Section [3] discusses the spectral prop- 
erties of our sample. In Section [5] we compare AMI/ VSA 
flux densities with WMAP flux densities. Finally, we sum- 
mar ize our results in Section HI] In a n acc ompanying pa- 
per, lAMI Consortium: Franzen et al.l l|2009l) (hereafter Pa- 
per II), we present results on the flux density variability of 
the sources. 



1.1 The NEWPS Catalogue 

The number of sources in the NEWPS catalogue was max- 
imised by exploiting low-frequency data. In the first in- 
stance, a catalogue was created from all those sources 
with S ^ 500 mjy in d ata at 4.85 GHz from the PMN 
{Griffith fc Wright! Il993l ) and GB6 {Gregory et al l Il996l ) 
surveys. Where there are holes in the sky coverage of these 



surv eys they were filled using data from the NVSS sur- 
vey {Condon et al.lll99Sl) at 1.4 GHz and the SUMSS sur- 
vey ( Bock. Large fc Sadler! [l999l ) at 843 MHz, again using 
sources with S ^ 500 mjy. The WMAP data were then fil- 
tered u sing the second member of th e Mexican Hat Wavelet 
family ( Gonza lez-Nuevo et alj [20061 ) at the positions of all 
the sources in the source catalogue. Any sources detected at 
^ 5a, in any one of the frequency channels, were included 
in the NEWPS 5-a catalogue. 

A small number of sources detected in the original 
WMAP analysis were not found in the input catalogue. This 
is because they have strongly inverted spectra, with flux 
densities < 500 mjy in the low-frequency data. The WMAP 
data were also filtered at the positions of these sources and, 
if detected at ^ 5a, were added to the catalogue. In total, 
at 33 GHz, 224 sources were detected at ^ 5a, five of which 
did not appear in the input catalogue. 



1.2 Source sample 

We selected sources in the NEWPS catalogue detected at 
> 5a at 33 GHz and with -4° sC 8 ^ 60° . These crite- 
ria were met by a total of 99 sources. It was not possible 
to obtain a good observation of J0528+2133, which hap- 
pens to lie just 1?5 from Tau A, with the VSA. We also 
excluded J2153+4716, which was identified as an HII region 
{Kuchar fc Clarkll 19971 ) . Having removed these two sources, 
we were left with a sample containing 97 sources and com- 
plete to 1.1 Jy at 33 GHz. 



2 THE TELESCOPES 



2.1 AMI 



The AMI (see AMI Consortium: 

lAMI Consortium: Zwart et all 120081 . for a full descrip- 
tion of the instrument) consists of two separate telescopes, 
the Large (AMI-LA) and Small (AMI-SA) Arrays; all the 
AMI data presented in this paper were obtained from the 
AMI-SA. The AMI-SA operates at frequencies between 
13.9 and 18.2 GHz with the passband divided into six 
channels of 0.72-GHz bandwidth. The average frequencies 
of the 6 usable bands are given in Table [T] (the AMI-SA 
has two frequency channels below 13.9 GHz which are not 
routinely used, owing to interference from satellites). The 
centre frequency is ~ 16.1 GHz. The results presented in 
this paper are a combination of continuum and channel flux 
densities. We make it clear in the results sections which of 
these we are using. 

The primary beam of the telescope at 16 GHz is ~ 20 ar- 
cmin FWHM. The synthesised beam, which is an effective 
measure of the resolution of the telescope, varies with fre- 
quency channel and observation declination. However, it is 
typically « 3 arcmin FWHM. The telescope measures a sin- 
gle, linear polarisation: Stokes parameter I + Q. It has a flux 
sensitivity of ~ 30 mjy in one second and is able to observe 
at declinations > —15°. 
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Table 1. Assumed flux densities for sources used for AMI flux 
density calibration. Note that the individual channel errors are 
not independent and that the error on each calibration error is 
about half a per cent. 
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2.2 VSA 

The main array of the VSA (see IWatson et al. I l2003l . for a 
detailed description of the instrument), which was used to 
obtain the VSA data presented in this paper, has a single fre- 
quency channel centred at 33 GHz with 1.5 GHz bandwidth. 
The VSA has several observing configurations. The observa- 
tions presented here were made with the VSA in its s uper- 
extended configuration (see iGenova-Santos et al1l200&l ). 

The flux sensitivity of the telescope in the super- 
extended configuration is ~ 2.7 Jy in one second. The pri- 
mary beam is 72 arcmin FWHM. Again, the synthesised 
beam varies with observation declination. A typical value, 
however, is 7 arcmin FWHM. The VSA can observe in the 
declination range —4° ^ 8 ^ 60°. The telescope is not 
equatorially mounted, so the linear polarisation measured 
changes with the hour angle of observation. At a parallactic 
angle of 0° the telescope measures Stokes parameter I — Q. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 AMI 

Observations of the WMAP sources were made using the 
AMI-SA during fourteen separate observing runs between 
2007 April and 2008 September. The majority of the sources 
were observed, at irregular intervals, at least three times dur- 
ing this period. During the earlier observing runs, which typ- 
ically lasted about 48 hours, the sources were observed twice 
(at different hour angles) for ten minutes, although some of 
the fainter sources were observed for longer to match the 
observing time on the VSA. This two-observation scheme 
provided a useful check of the reliability of our calibration 
procedures (see below). The two datasets were then con- 
catenated so as to improve the signal-to-noise and the uv 
coverage of the aperture plane. Having established the re- 
liability of our calibration procedures, during the later ob- 
serving runs the sources were observed only once, typically 
for twenty minutes each. 

Data reduction was carried out using reduce, our soft- 
ware developed, originally, for the VSA and later modi- 
fied for the AMI. This applies path delay corrections, auto- 
matic flags for intereference, pointing errors, shadowing and 
hardware faults, applies phase and amplitude calibrations, 
Fourier transforms the data into the frequency domain, and 



writes it out in uv fits format for imaging in the AlPij]] 
package. 

Flux density calibration was performed using short ob- 
servations of 3C48 and 3C286 interspersed with the obser- 
vations of the WMAP point sources. The flux densities as- 
sumed for these sources, used for calibrating each of the 
freque ncy channels, are consistent with those o f B aars et alj 
l| 19771 ). Sinc e we measure a different polarisation from that 
measur 

edbyH aars et all (I + Q, as opposed to I), we correct 
for the polarisations of the calibrator sources by interpolat- 
ing VLA data collected at 5, 8 and 22 GHz. The assumed 
flux densities in each of the AMI frequency channels are 
listed in Table [1] 

We correct for changing air mass and variations in at- 
mospheric conditions both during the observations and be- 
tween the observations of the sources and the flux density 
calibrators by monitoring a modulated noise signal, injected 
at the front-end of each antenna. The data are also weighted 
by comparing the ratio of the power of the modulated noise 
signal to the total power input to the correlator (which is 
kept constant), to that obtained in cool, dry, clear weather 
conditions. Samples are flagged if the data were taken in 
poor weather conditions, or if there were large changes in 
the weather between the observation of the source and the 
flux density calibrator. 

The telescope is not expected to be phase stable over 
the duration of an entire observing run. It is, however, phase 
stable over the length of an individual observation and, since 
the SNRs for all sources were sufficiently large, we were able 
to use them as their own phase calibrators for the purpose of 
absolute phase calibration. This negated the requirement for 
observations of interleaved phase calibrators. We find that 
implementing this scheme for phase calibration produces a 
significant decrease in the dispersion of the flux densities 
measured for the cross-calibrated observations of 3C286 and 
3C48. It also reduces the scatter in the two flux densities, 
measured for each of the WMAP point sources, within the 
same observing run. 

The uv fits data output from reduce were imaged 
using the AIPS package. Maps were produced for channels 3 
to 8 individually and combined. The great majority of the 
sources were not resolved by the AMI-SA. In these cases, the 
quoted flux density is the peak value from the dirty map. 
For the few sources that were resolved we first CLEAN the 
map and then measu re the integr ated flux density. Following 
a similar method to iReesI i|l990h . we sum contiguous pixels 
down to a contour level of half the peak flux density value 
to give Imp- We then apply the same operation to the beam 
with a contour level of half the height of the maximum on 
the beam, to give hm- The integrated flux density is then 
taken to be I mp /Ibm.- 

3.1.1 Flux density error estimation 

We checked the flux density calibration of the telescope by 
comparing the two observations of each source from individ- 
ual observing runs, making the assumption that the vari- 
ability timescale was long compared to the duration of the 
observing run. We also cross-calibrated observations of our 

1 http://www.aips.nrao.edu 
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flux density calibrators. Both tests indicated that the flux 
density calibration of the telescope is consistent to « 3 per 
cent r.m.s. 

We have taken a conservative approach and assumed 
that the error, a, on a measured continuum flux density, S, 
consists of a 4 per cent calibration error added in quadrature 
with the thermal noise on the map, n (which is measured 
far from the primary beam). We add the calibration error in 
quadrature with the thermal noise because these two errors 
are uncorrelated. The total error is therefore given by a — 
\J (0.04S) 2 + n 2 . Similar tests were performed to measure 
the calibration errors on the individual frequency channels, 
which are shown in Table [T] Errors quoted on channel flux 
densities consist of the map noise added in quadrature with 
the calibration error for the relevant channel. 



3.2 VSA 

Observations of the WMAP sources were performed dur- 
ing ten separate observing runs between 2007 February and 
2008 July. The majority of the sources were observed, at 
irregular intervals, at least four times during this period. 
Observations were 10-20 minutes in length. During most 
observing runs, which typically spanned 24-48 hours, each 
of the sources was observed twice. A small proportion of 
observations was discarded for reasons unrelated to source 
characteristics, such as interference from the Sun, Moon or 
planets. 

A full description of t he VS A calibration process is pre- 
sented in IWatson et al.l l|2003h . Briefly, the absolute flux 
density calibration of VSA data is determined from ob- 
servations of Jupiter, the brightness temperature (Tj up = 
146.6 ± 0.7 K) of which is taken from the WMAP 5-year 
data. These, in turn, were calibrated on the CMB dipole 
l|Hill et al.l 120081 ') . This flux density scale is transferred to 
our other calibrator sources: Tau A, Cas A, Cyg A and Sat- 
urn. The specifications for the super-extended configuration, 
such as the correction for the fact that Tau A and Cas A are 
partially resolved in the longest baselines, are essentially the 
same as those a dopted for the e x tende d configuration and 
are explained in iDickinson et all (|2004h . The assumed flux 
density of Cyg A, which is unresolved by the VSA, is 36.7 Jy 
at 33 GHz. 

Data reduction was carried out using the reduce pack- 
age written specifically for the VSA (see Section [3.1J1 . Am- 
plitude and phase calibrations were performed using obser- 
vations of Cas A, Cyg A, Saturn and Tau A interspersed 
with the observations of the WMAP sources. A calibra- 
tion scheme addres sing tempe rature effects on the phase was 
adopted (see lLancasterl 120041 ) . The data were corrected for 
changes in system temperature and in atmospheric opacity 
with elevation, based on the monitoring of modulated noise 
signals injected into the VSA system. The correction was 
typically a few per cent. In order to achieve the optimum 
overall noise level, the data were re-weighted based on the 
r.m.s. noise of each baseline. 



3.2.1 Phase calibration 

The uv fits data output from reduce were imaged using 
the MPS package. To correct phase errors, we self-calibrated 
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Figure 1. The distribution of the ratios, p, of the flux densities 
with self-calibration to those without for the 176 observations 
with sufficient SNR for the CALIB task to be applied successfully. 



datasets making the assumption that all the sources are 
unresolved by the VSA. To achieve this, we employed the 
CALIB task using a point-source model. We used the phase- 
only solution mode, with a solution interval encompassing 
the entire observation. We ran the task on each dataset but 
only applied solutions to the data if a solution was found 
for each antenna. We found this to be the case for a total of 
176 observations (~ 25 per cent of observations) , which typ- 
ically had signal-to-noise ratios (SNRs) > 15. The remaining 
observations had insufficient SNR for self-calibration to be 
applied successfully (we explain how we addressed this prob- 
lem in the next paragraph). We then produced dirty maps 
using the IMAGR task, applying natural weighting in order to 
maximise the SNR. In each map, we extracted the peak flux 
density inside a square at the map centre with half- width of 
10 arcmin. 

For the 176 observations with sufficient SNR for the 
CALIB task to be applied successfully, we also produced maps 
using the un-self-calibrated data in order to compare flux 
densities with and without self-calibration. Fig. [TJ shows 
the distribution of the ratios, p, of the flux densities with 
self-calibration to those without. We emphasize that we are 
dealing with peak flux densities. As expected, phase errors 
systematically result in flux densities being underestimated. 
The mean and median values of p are 1.048 and 1.036 re- 
spectively. In response to these findings, for the remaining, 
un-self-calibrated observations, we, accordingly, have multi- 
plied flux densities by 1.048 as a correction. 



3.2.2 Flux density error estimation 

The thermal noise was measured from the time-series dis- 
persion of the visibilities in reduce and is, for a typical 
observation, 0.15-0.20 Jy. We checked the flux density cali- 
bration of the telescope by cross-calibrating observations of 
our flux density calibrators. This test indicated that, from 
the amplitude alone, the flux density calibration of the tele- 
scope is subject to an error of « 4 per cent r.m.s. From 
the standard deviation of p (0.041), we deduced that, from 
the phase alone, the flux density calibration is subject to an 
error of ~ 4 per cent r.m.s. The fact that this error has a 
non-Gaussian distribution, as is clearly apparent in Fig. [1] 
should not be a significant problem, because, in the cases 
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where it is not possible to apply self-calibration, the total 
error is dominated by the thermal noise. 

We have, therefore, assumed that, for each source suf- 
ficiently bright to be used as its own phase calibrator, 
the total error on a measured flux density, S, is given by 
a = 1V (0.04S) 2 + n 2 , where n is the thermal noise. Other- 
wise, taking a conservative approach, the total error is taken 
to be o = \J (O-OSS*) 2 + n 2 . We add the calibration error in 
quadrature with the thermal noise because these two errors 
are uncorrelated. 

We checked our error estimates by comparing the two 
observations of each source from individual observing runs, 
making the assumption that the variability timescale was 
long compared to the duration of the observing run. From 
a total of 687 observations, there were 248 such pairs of 
observations. A x 2 -test for the difference between the two 
flux densities produced a reduced x 2 -value of 0.98 for 247 
degrees of freedom. The probability of exceeding this value, 
by chance, is 0.58. As a result, we are confident in our error 
estimates. 

For each source observed twice during the same observ- 
ing run, we have treated the two observations as a single 
observation. If self-calibration was applied to both observa- 
tions, the two datasets were concatenated before imaging 
and the quoted flux density, S, is the peak flux density in 
the resulting map. Since, for a self-calibrated source, the 
phase errors have been corrected, the quoted error is given 
by a = ^(O.OAS) 2 + n 2 , where n = (l/n\ + l/n|) ~ 1/2 , 
and tia and ub are the thermal noises in the two individual 
observations. 

Because of possible discrepancies in source positions 
arising from phase errors, we concatenated the two datasets 
only in cases where it was possible to apply self-calibration 
to both observations. Otherwise, the quoted flux density is 
given by 

_ Sa/u 2 a + S B /n 2 B 

1/ni + l/ni ' () 

where Sa and Sb are the two measured flux densities with 
thermal noises ua and n_g respectively. 

If self-calibration was applied to neither of the two ob- 
servations, the quoted error is given by a — yj (0.08S 1 ) 2 + n 2 
and in, the few cases where it was applied to one of the 
two observations, a = yj (0.06S) 2 + n 2 . Note that we have 
considered the worst case scenario where the flux density 
calibration errors in the two individual observations are cor- 
related (which is a possibility given how the observations 
were scheduled). 

3.2.3 Peak- flux- density bias 

We have identified a bias in our measurements at low SNR, 
which results from measuring peak flux densities as opposed 
to those at the exact source positions. Low-frequency iden- 
tifications are given for each of the sources in the NEWPS 
catalogue. For increased accuracy, we used the coordinates 
of sources listed as likely matches in the 4.85-GHz PMN or 
GB6 catalogues. Because of the VSA's very different observ- 
ing frequency, we checked for positional differences (these 
were typically 1.5 arcmin), but, in fact, these were mainly 
caused by phase calibration errors in the VSA data. There- 
fore, in each case, rather than measuring the flux density at 
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Figure 2. Results of a simulation to investigate the bias resulting 
from measuring peak flux densities at low SNR with the VSA. We 
simulated pointed observations of a 1.0-Jy point source with dif- 
ferent SNRs, performing 500 realisations per SNR. Open squares 
show the average flux densities measured at the pixel correspond- 
ing to the map centre; filled circles show the average peak flux 
density inside a square, centred at the map centre, with half-width 
of 10 arcmin. 

the pixel corresponding to the pointing centre, we measured 
the peak flux density inside a square at the map centre with 
half-width of 10 arcmin (i.e. within « one synthesised beam 
from the pointing centre). However, at low SNR, the peak 
flux densities are biased to be slightly high, because of the 
increased likelihood of peak positions lying on top of positive 
noise fluctuations. 

We have perfomed simulations to investigate this bias, 
the results of which are shown in Fig. [2] Pointed observa- 
tions of a 1.0- Jy point source with SNRs ranging between 
2.0 and 7.0 were simulated. We used a typical uv coverage 
for the VSA and produced dirty maps using the AIPS pack- 
age. We performed a total of 500 realisations per SNR. For 
each simulated map, we measured the flux density at the 
pixel corresponding to the map centre and the peak flux 
density inside a square at the map centre with half- width of 
10 arcmin. 

As expected, we found that there is no bias if flux densi- 
ties are measured at the correct positions, no matter how low 
the SNR. We found that, at SNR = 5.0, peak flux densities 
are on average « 4 per cent higher than the true flux den- 
sity. The effect becomes more severe as the SNR decreases 
further, but even at SNR = 3.0 (only as 3 per cent of ob- 
servations are below this SNR), peak flux densities are on 
average only « 12 per cent higher than the true flux densi- 
ties. Since the magnitude of the bias is small relative to the 
thermal noise, we do not expect it to have any significant 
effect on eventual results and have, therefore, not corrected 
the VSA flux densities for this effect. 



4 13.9 TO 33.75 GHZ SPECTRAL PROPERTIES 

In Fig. [3]we have plotted the mean flux density at 16 GHz 
versus that at 33 GHz for each of the sources in our sample. 
The plot demonstrates that the sample contains a large pro- 
portion of sources with rising spectra. The mean AMI and 
VSA flux densties, along with the number of observations, 
are given for each source, in Table [2] 
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Figure 3. Mean flux density at 16 versus that at 33 GHz. The 
dashed line corresponds to = 
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Figure 4. Histogram of the spectral index distribution o?2 I 5 . 



We have simultaneous observations of most of the 
sources in the sample, a significant advantage over other 
work. The obervations on the two telescopes were typically 
performed within a few minutes of one another; at worst 
they were separated by two days. We have at least a sin- 
gle pair of reliable, simultaneous observations for 94 of the 
sources. For the remaining three sources in our sample we 
have no reliable simultaneous observations; this is due to 
hardware problems or interference during the observations. 

We have fitted power-law spectra to these simultaneous 
data, defining spectral index, a, by S oc v~ a . In fitting the 
spectra we have made use of the spectral data from the six 
useable AMI frequency channels and the single VSA chan- 
nel. The frequencies at the lower end of the AMI and the 
upper end of the VSA passbands are 13.9 and 33.75 GHz 
respectively. We, therefore, quote spectra, afjjig 5 , between 
these two frequencies. 

The spectra were fitted by sampling using a Markov 
Chain Monte Carlo technique from a uniform likelihood 
function. This method takes into account the asymmetric 
errors in In S and enables the calculation of an error esti- 
mate for each spectral index, taking account of the posterior 
distribution. We find that, if there were to be a 5 per cent 
difference in the flux density scales between the two tele- 
scopes, the resulting systematic error in the spectral index 
would be 0.07. 

We provide a histogram of the spectral index distribu- 
tion aJllg 5 for the 94 sources in Fig. [4] For those sources 
with multiple pairs of simultaneous observations, we have 
used the mean spectral index in plotting the histogram. The 
spectral index, fitted to the simultaneous observations, of 
each source, along with the associated error, are listed in 
Table [21 The table shows that the typical errors are smaller 
than the bin size used in plotting the histogram. In Fig. [5]we 
have provided spectra for a sample of typical sources. Of the 
94 sources, 41 (44 per cent) have rising spectra (afjjig 5 < 0) 
and 87 (93 per cent) have spectra with ailig 5 < 0.5. The 
median spectral index is 0.04. 



5 COMPARISON OF AMI/ VSA AND WMAP 
FLUX DENSITIES 

In Fig. [6] we compare our 33-GHz flux densities with the 33- 
GHz flux densities measured by WMAP and, in figure [7] we 
compare our 16-GHz continuum flux densities with the 23- 
GHz flux densities measured by WMAP. Because the sources 
are variable (see Paper II), we used the mean AMI/ VSA 
flux density for each source. Since the level of variability for 
many of the sources is large compared to the errors on the 
flux density measurements, we have not included error bars 
in these figures. 

We have use d the Bayes-corrected version (see 
iHerranz et"afll2006l ) of the NEWPS Catalogue to take into 
account Eddington bias. Close to the survey completeness 
limit, source flux densities are biased high. The presence of 
noise, and the fact that faint sources are more numerous 
than bright ones, result in an increased likelihood for peak 
positions to lie on top of positive fluctuations. The correc- 
tion for each source in the NEWPS catalogue is calculated 
from the SNR and the slope of the source counts. The source 
counts are described by a power law, which is estimated from 
the NEWPS flux densities themselves. 

The flux density calibrations of all the WMAP chan- 
nels have been reassessed through estimates of the effec- 
tive beam areas and correct ion factors have been derived 
l|Gonzalez-Nuevo et alj 12008). These correction factors are 
1.05, 1.086, 1.136 and 1.15 at 23, 33, 41 and 61 GHz re- 
spectively. We have corrected the NEWPS flux densities by 
these self-calibration factors. 

We note that eight sources have no flux densities quoted 
at 23 GHz in the NEWPS 5-cr catalogue because they are 
not detected at 5a. For these sources, we used the flux 
densities quoted in the NEWPS 3-cr catalogue. These sources 
are represented as open squares in Fig. [7] 

Figures. [6] and [7] indicate that there is generally very 
good agreement between the AMI/ VSA flux densities and 
the WMAP ones. Discrepancies between the flux densities 
are expected given our findings in Paper II concerning the 
general level of variability of the sources in the sample. Our 
observations were carried out in 2007 and 2008, i.e. in a pe- 
riod after that of WMAP 3-year maps (obtained by averag- 
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Figure 5. Examples of spectra between 13.9 and 33.75 GHz. Note that the errors on the AMI flux densities measured in different 
frequency channels are not independent. 



ing over the data collected during 2001-2004). We note that 
the frequencies at which the AMI and WMAP results are 
compared are not the same. However, we do not expect this 
to contribute significantly to the scatter in Fig. [7] because 
the two frequencies are relatively close to one another and, 
as discussed in Section [4] a large proportion of the sources 
have flat spectra. 

We have checked whether there is a good agreement be- 
tween the VSA and WMAP flux densities for a very bright 
non-variable source in our sample, Cyg A (J 1959+4034). 
This source is also one of the VSA calibrator sources, 
with an assumed flux density of 36.7 Jy. The overall ex- 
tent of the source is 2.1 arcmin and it is unresolved by 
the VSA and WMAP. No variability in the core of Cyg A 
has been reported, at any frequency. In any event, the 
core represents a small fraction of the total flux density. 
lAlexander. Brown fc Scottl (|l984l 'l made observations of Cyg 
A with the Cambridge 5-km radio telescope and measured 
the flux density of the central component at 15.4 GHz to be 
just 1.22 ±0.2 Jy. We found that the flux density quoted by 
WMAP (37.87 ± 0.49 Jy) is in excellent agreement with the 
flux density measured, here, by the VSA (37.43 ± 0.84 Jy). 

5.0.1 Offset between AMI/VSA and WMAP flux densities 

It is also apparent in figs. [6] and [7] that there is a system- 
atic offset at faint flux densities, the AMI/VSA flux densi- 
ties tending to be lower than the WMAP ones. We quan- 
tified this effect as follows: for each source, we calculated 



the mean AMI/VSA and WMAP flux density and on this 
basis, split the sources into three equally sized groups, a 
low-, medium- and high-flux-density group. In each group, 
we then compared the mean AMI/VSA flux density with the 
mean WMAP flux density. Given the large range of flux den- 
sities, we performed calculations in log space. The results are 
shown in Figs. [8] and [9] The error bars are standard errors 
of the means. 

In the high-flux-density group, the mean AMI/VSA and 
WMAP flux densities are in excellent agreement. However, 
in the low-flux-density group, the mean AMI and VSA flux 
densities are significantly lower than the mean WMAP flux 
density, the discrepancies being significant at a level of 3.0a 
(AMI) and 4.9cr (VSA). In the medium-flux-density group, 
the effects are significant at levels of 4.2a (AMI) and 2.0a 
(VSA). This statistical analysis suggests that the systematic 
offset at faint flux densities is very unlikely to be due to 
chance and requires an explanation. 

One possible issue is that the AMI beam is « 3 arcmin, 
the VSA beam » 7 arcmin, the WMAP beam in the K- 
band « 53 arcmin and the WMAP beam in the Ka-band 
w 40 arcmin. However, incomplete sampling of extended 
sources with AMI and the VSA cannot explain the bias be- 
cause, as discussed in Section[37l] practically all sources were 
found to be point-like at the AMI resolution. 

Another possible factor is source confusion with 
WMAP. In order to have a significant effect, this would re- 
quire more than one source with S > 1 Jy to lie within a 
WMAP resolution element. This is very unlikely given the 
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Figure 6. Comparison of VSA 33-GHz flux densities with 
WMAP 3-year 33-GHz flux densities, with a line indicating equal 
flux density values. 
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Figure 7. Comparison of AMI 16-GHz flux densities with WMAP 
3-year 23-GHz flux densities, with a line indicating equal flux den- 
sity values. We used flux densities quoted in the NEWPS 3<r cat- 
alogue for the faintest sources that did not appear in the NEWPS 
5-cr catalogue. These sources are represented as open squares. 



number of sources with S > 1 Jy detected by WMAP in 
all-sky maps. The NEWPS catalogue contains 224 sources 
at 33 GHz; there are ~ 500 WMAP beam areas per source 
in the Ka-band. 

Given the length of time between our observations and 
the WMAP ones and the general level of variability (see 
Paper II), we suggest that the overwhelming cause of the 
systematic effect is Eddington bias arising from variability. 



5.0.2 Eddington bias in WMAP 5-year data 

We have used the WMAP 5-year, 33-GHz data to check 
whether the above effect is present in the WMAP data. We 
ran a blind search with a 4a cut on the first-year map in the 
region 20° ^ b 90° ; we excluded sources in the region 0° ^ 
b < 20° to avoid galactic contamination. We followed this by 
a non-blind search on the five single-year maps, using prior 
knowledge of source positions from the first-year map. We 
were then left with a list of sources selected in the first-year 
map, with flux densities in each separate year. We corrected 
flux densities from the first-year map for the Eddington bias. 
The correction for this bi as becomes unreliable at low SNRs 
(see iHerranz et all [2006). For this sample of sources, the 
correction starts becoming unreliable at SNR < 4.1. Having 
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Figure 8. Comparison of the mean VSA 33-GHz flux density 
with the mean WMAP 3-year 33-GHz flux density in three groups 
of sources, a low, medium and high flux density group, with a line 
indicating equal flux density values. The error bars are standard 
errors of the means. Calculations were performed in log space. 
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Figure 9. Comparison of the mean AMI 16-GHz flux density 
with the mean textitWMAP 3-year 23-GHz flux density in three 
groups of sources, a low, medium and high flux density group, 
with a line indicating equal flux density values. The error bars 
are standard errors of the means. Calculations were performed in 
log space. 



rerun the source-finding on the first-year map in the non- 
blind approach, some of the first-year flux densities were 
found to be significantly below 4cr. We, therefore, removed 
sources detected below 4.1cr in the first-year map. 

In Fig. 1101 we compare flux densities in year 1 with mean 
flux densities in years 2 to 5. We followed the same statistical 
analysis as described above, the results of which are shown in 
fig. 1111 As expected, in the low-flux-density group, the mean 
flux density in year 1, in which the sources were selected, is 
significantly higher, the effect being significant at a level of 
3.5 a. There is no apparent bias in the medium- and high- 
flux-density groups. 

In order to verify that these results are not simply due to 
changes in telescope systematics over the course of the 5-year 
survey, we repeated the analysis selecting sources in the fifth- 
year map instead. This time, we selected sources from the 
whole sky. We removed sources detected below 4.1rj in the 
fifth-year map as well as sources in the strip \b\ ^ 20°. Fig. 1121 
compares flux densities in year 5 with mean flux densities in 
years 1 to 4 and Fig. [13] shows the results of our statistical 
analysis. We note that the error bars in Fig. [11] are bigger 
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Figure 10. Comparison of flux densities in year 1 with flux densi- 
ties in years 2 to 5 for WMAP 33-GHz data. Sources were selected 
in the first-year map, with a line indicating equal flux density val- 
ues. See text for more details on the analysis. 
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Figure 12. Comparison of flux densities in year 5 with flux densi- 
ties in years 1 to 4 for WMAP 33-GHz data. Sources were selected 
in the fifth-year map, with a line indicating equal flux density val- 
ues. See text for more details on the analysis. 
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Figure 11. Comparison of the mean WMAP 33-GHz flux density 
in years 2 to 5 with the mean WMAP 33-GHz flux density in 
year 1 in three groups of sources, a low-, medium- and high-flux- 
density group, with a line indicating equal flux density values. 
The error bars are standard errors of the means. Calculations 
were performed in log space. 
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Figure 13. Comparison of the mean WMAP 33-GHz flux densi- 
ties in years 1 to 4 with the mean WMAP 33-GHz flux densities in 
year 5 in three groups of sources, a low-, medium- and high-flux- 
density group, with a line indicating equal flux density values. 
The error bars are standard errors of the means. Calculations 
were performed in log space. 



than those in 1131 This is simply because the analysis was 
carried out using data obtained from sky areas of different 
size. 

In the low-flux-density group, the mean flux density in 
year 5, in which the sources were selected, is significantly 
higher, the effect being significant at a level of 6.9 a. Again, 
there is no apparent bias in the medium- and high-flux- 
density groups. We conclude that the bias is clearly present 
in the WMAP 5-year data and that it is not related to the 
noise or to telescope systematics, but rather is a consequence 
of source variability. In any survey, there will be preferen- 
tial selection of those sources which were above their mean 
flux density values at the time of observation rather than 
below. The magnitude of the bias will not just depend on 
the general level of variability and the slope of the source 
counts, but also the survey duration and the typical vari- 
ability timescale. It should gradually decrease as the survey 
duration (i.e. making repeated measurements of the source 
flux density over a period of time) becomes comparable to 
the typical variability timescale. 



6 CONCLUSIONS 

In order to tie down the cm-wave spectra of sources at 
the high-flux-density end of the source population at cm 
wavelengths, we have made observations with the AMI- 
SA (13.9-18.2 GHz) and the VSA (33 GHz) of a complete 
sample of sources found with WMAP at 33 GHz and with 
S33 GHz ^1.1 Jy; AMI and VSA observations were scheduled 
so that variability is not an issue in the spectral measure- 
ments. We found that: 

(1) the spectra are very different from those of bright 
sources at low frequency: 93 per cent have spectra with 
«?l:9 5 < 0.5; 

(2) 44 per cent have rising spectra (ailig 5 < 0.0); 

(3) in the group of 33 sources with average WMAP and 
VSA flux densities > 2.26 Jy, the average flux densities mea- 
sured by WMAP and VSA are fully consistent with each 
other. However, in the group of 32 sources with average 
WMAP and VSA flux densities < 1.36 Jy, the mean VSA 
flux density is lower than the mean WMAP flux density, the 
discrepancy between the two being significant at a level of 
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4.9a. We ascribe this to Eddington bias arising from vari- 
ability. 
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Table 2. Results for individual sources. The columns are the source name, taken from 
the NEWPS catalogue; the Equatorial coordinates (J2000), from the PMN or GB6 
catalogues; the number of observations, n; the mean flux density, S; and the spectral 
index, ofg-g 5 . 
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Table 2 — continued 
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